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Abstract 

The observation of a sharp spectral feature in the gamma-ray sky would be one of the cleanest 
ways to identify dark matter and pinpoint its properties. Over the years a lot of attention has 
been paid to two specific features, namely gamma-ray lines and internal bremsstrahlung. Here, we 
explore a third class of spectral signatures, box-shaped gamma-ray spectra, that naturally arise 
in dark matter cascade annihilations or decays into intermediate particles that in turn decay into 
photons. Using Fermi-LAT data, we derive constraints on the dark matter parameter space for 
both annihilating and decaying dark matter, and show explicitly that our limits are competitive to 
strategies employing standard spectral features. More importantly, we find robust limits even in 
the case of non-degenerate dark matter and intermediate particle masses. This result is particularly 
relevant in constraining dark matter frameworks with gamma-ray data. We conclude by illustrating 
the power of box-shaped gamma-ray constraints on concrete particle physics scenarios. 



I. INTRODUCTION 



The quest for dark matter (DM) is entering a crucial period as collider, direct and indirect 
searches close in on the most popular frameworks for weakly interacting massive particles 
(WIMPs). No clear, unambiguous signal has been pinpointed as yet and correspondingly 
stringent constraints on the properties of dark matter have been derived. Indirect searches, 
in particular, offer a multitude of channels and potential targets as discussed at length in 
the literature (see for instance [IH1] for reviews). However, indirect signals are, on the one 
hand, prone to different sorts of uncertainties (including cosmic-ray propagation [oTfTU] and 
the dark matter distribution in our Galaxy [TTVTH] ) and, on the other hand, lack strong 
discrimination power against astrophysical sources [T5(, [16] . This makes any claim of iden- 
tification of dark matter through indirect channels - as well as the reconstruction of the 
underlying properties - extremely problematic. In this context, gamma-ray spectral fea- 
tures assume today a particularly relevant role since there is no known astrophysical source 
capable of mimicking such signature. Moreover, the LAT instrument [T7J HH] on board the 
Fermi satellite has dramatically improved our mapping of the gamma-ray sky and offers 
unprecedented sensitivity in the search for spectral features. Therefore, gamma-ray signa- 
tures rest as perhaps the best indirect strategy to identify dark matter in an unambiguous 
fashion. Signatures such as gamma-ray lines at (a fraction of) the dark matter particle 
mass rriDM an d the hard spectra provided by internal bremsstrahlung just below m^M have 
been thoroughly investigated [IHH2I] and shown to place strong limits on particle physics 
models (for a tentative gamma-ray line, see the recent works j25j EE]). Here, we focus on 
a different spectral feature that is somewhat unexplored in the existing literature. If dark 
matter annihilates or decays into an intermediate particle <fi that in turn decays appreciably 
into photons, the resulting photon spectrum resembles a "box" whose centre and width are 
entirely fixed by m DM and m^. The aim of the present paper is to study the relevance of 
these box-shaped features for dark matter searches and to show that they can be at least as 
constraining as the widely-studied gamma-ray lines. We start by discussing the phenomenol- 
ogy of dark matter cascade processes and the corresponding photon fluxes in Section [TT| In 



Section III we derive our constraints on the dark matter parameter space using Fermi-LAT 
data and briefly address the detectability of box-shaped gamma-ray features. Finally, we 
present concrete particle physics models for cascade annihilations into photons in Section 
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IV 



before concluding in Section |V| 



II. GAMMA-RAY SPECTRAL FEATURES FROM CASCADE ANNIHILATIONS 
OR DECAYS 

Dark matter cascade annihilations into light degrees of freedom have been explored before 
[27H22] as a means of yielding sizable fluxes of electrons and positrons, diffuse gamma- 
rays, synchrotron radiation and neutrinos. We focus instead on the possibility of producing 
distinctive gamma-ray spectral features through 1-step cascade processes. Suppose that dark 
matter self-annihilates into a pair of scalars that in turn decay into a pair of photons 1 



(for concrete particle physics realisations and branching ratios see Section IV). Each of the 
four photons emitted per annihilation has a monochromatic energy E' = in the rest 

frame of the corresponding scalar 0. In the lab frame - where the dark matter particles are 
non-relativistic and the scalars have energy E& = ttidm - the photon energy reads 



lm DM \ y m DM 

with 9 the angle between the outgoing photon and the parent scalar in the lab frame. From 
this equation we read that the spectrum has sharp ends defined by the parameters m^M 
and m<f). The highest (lowest) energy corresponds to a photon emitted at an angle 8 = 0° 
(180°) with respect to the momentum of the parent scalar. Since the decaying particle is a 
scalar, the photon emission is isotropic. Hence the resulting spectrum is constant between 
the energy endpoints and takes a flat, box-shaped form: 

^.^J-e {E - E _ M E + -E) , (2) 



where 6 is the Heaviside function, AE = E + — E_ = J m 2 DM — m| is the box width and 
E± = (m DM /2) ^1 ± — ra 2 ^/ra 2 DA/ ^j are the box edges. In the case of dark matter decays 
into a pair of scalars, the above expressions apply with the replacement moM ~^ m DM /2. 



1 The case of dark matter self-annihilation into one scalar and another particle and/or the decay of the 
scalar into one photon and another particle is also feasible, but leads to very similar phenomenology as 
presented here and is therefore omitted. 
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Let us briefly discuss the impact of the primary spectrum ^ on gamma-ray searches 
for dark matter. First of all, notice that the centre of the "box" E c = (E + + EJ)/2 
lies at mDAf/2, while the relative width is chiefly regulated by the degeneracy parameter 



Am/m DM = {rriDM - rn<{,)/m DM G [0,1] since AE/E C = 2 v / 2Am/m DM - Am 2 /m 2 DM . 
Hence, ttidm and Am/mDM fix the shape of the boxy spectrum, and accordingly we shall 
use both as our main phenomenological parameters. Secondly, in the limit Am/mDM 
(or, effectively, when AE/E C falls bellow the experimental resolution) the spectrum dN 1 /dE 1 
reduces to a line. But, unlike the well-known 77 and 7Z lines, in this case there are four 
photons emitted per annihilation and each carries half of the dark matter particle rest-mass 
energy, or one-fourth for decaying dark matter models. Clearly, as we consider less and less 
degenerate configurations, i.e. as Am/mDM 1, the gamma-ray box gets wider in energy 
and dimmer in amplitude which might lead one to conclude that non-degenerate cases are 
hardly observable. However, as we shall see, even for Am/mDM — 1 the spectral plateau sits 
at non-negligible amplitudes and its feebleness is partially compensated by the extension to 
high energies. In fact, given that all known gamma-ray backgrounds (and also some spectral 
features such as internal bremsstrahlung) fall rapidly with energy, a flat spectrum at high 
energies - albeit dim - has good detectability prospects and strong constraining power. This 



is the backbone of our approach and shall be explicitly illustrated in Section III 



Now, the flux of photons at the Earth induced by the annihilation spectrum ^ in the 
case of Majorana dark matter particles is given by [TJ [21 [23] 

07( 7) = dE~dSdmt ~ dE~, M J An ann ' { } 

where (av) is the thermally averaged annihilation cross-section into <fi pairs, AQ is the 
observed field of view and J ann = j, ds p 2 DM is the integral of the squared dark matter 
density pdm along the line of sight. For models where the dark matter particle is not 
Majorana, equation ^ is to be multiplied by a factor 1/2. The photon flux induced by 
decaying dark matter is obtained by replacing (av ) /8iim DM — > Y/AnmoM (r is the decay 
rate) and J ann — > Jd ec = f los ds pdm in equation (|3]), as well as in dm ~ > wiDjy/2 in equations 
g) and @. In any spectral analysis it is essential to convolute the expected flux 7 with 
the experimental energy resolution. Assuming a gaussian detector response characterised 
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by a standard deviation cr(E), the convoluted photon flux reads 



dE' 4>,(E'[ 



exp 



(4) 



This procedure smears all sharp edges and spikes eventually present in 7 and ultimately 
limits the search for spectral features. 
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FIG. 1: The spectrum of box-shaped gamma-ray features. The left panel displays the unconvoluted 
(dashed) and convoluted (solid) box spectra for uidm = 100 GeV, = 60 GeV and (av) = 3 x 
10 -26 cm 3 /s on top of the 2-yr Fermi-LAT data (borrowed from |23| ) for the galactic centre region. 
The right frame shows the convoluted box spectra for thdm = 100 GeV, (av) = 3 x 10~ 26 cm 3 /s 
and several values of m^. 



We are interested in exploring box-shaped gamma-ray features in the energy range rele- 
vant for typical WIMPs, i.e. from a few GeV to a few TeV, so we shall focus on Fermi-LAT 
performance and data all through the manuscript. The energy resolution of the LAT instru- 
ment is assumed gaussian and parameterised according to [33, El], giving a(E)/E = 8 (12)% 
at E = 1 (200) GeV. Considering a more realistic non-gaussian LAT energy response would 
lead to a slightly different flux drop at the edges of the box spectrum; that would be rele- 
vant for the identification of a spectral feature but is of reduced importance in deriving flux 
constraints as we pursue in the present article. We consider as our main regions of interest 
the galactic centre and halo regions as defined in [23] (cf. Table 1 therein). The former 
features Aft = 1.30 sr, J AQ dVt J ann = 9.2 x 10 22 GeV 2 cm" 5 sr and J An dVt J dec = 6.9 x 10 22 

5 



GeV cm -2 sr, while the latter presents Af2 = 10.4 sr, J An d^lJ an n = 8.3 x 10 22 GeV 2 cm -5 
sr and J An dVL J dec = 2.2 x 10 23 GeV cm -2 sr, assuming a Navarro-Frenk-White (NFW) 
profile normalised to a local dark matter density of 0.4 GeV/cm 3 . Following the findings 
of [23], we shall focus on the centre (halo) region to derive constraints on annihilating (de- 
caying) dark matter. For the centre region, figure [I] (left) shows the unconvoluted and 
convoluted box spectra taking uidm = 100 GeV, = 60 GeV (or Am/rriDM = 0.4) and 
(av) = 3 x 10~ 26 cm 3 /s, as well as the 2-yr Fermi-LAT data borrowed from the analysis 
in [23] and distributed into 20 energy bins per decade. We have checked a posteriori that 
our results do not vary significantly with the particular choice of the binning. The reader 
is referred to [23] for further details on the Fermi-LAT data used here. Figure [I] (right) il- 
lustrates instead the effect of varying the mass degeneracy parameter Am/rriDM- The plots 
highlight the key phenomenological features of the dark matter models under scrutiny. As 
discussed above, in the limit of vanishing Am the box spectrum squeezes to a gamma-ray 
line whose amplitude and width are fixed by the energy resolution cr(E). Furthermore, the 
right panel of figure [T] clearly shows the relative importance of a flat spectrum extending to 
high energies. Notice as well that the convoluted box is not exactly symmetric about its 
centre because of the energy-dependent a(E)/E. Since the gamma-ray flux measured by 
Fermi-LAT is softer than E~ 2 , we can already anticipate that the most constraining energy 
bin for a given box spectrum is usually the one at which -E 2 7 peaks. 

Before proceeding with the derivation of Fermi-LAT constraints, a few comments are 
in order. The shape of the dark matter profile entering in the J— factors is not precisely 
known as of today despite the ground-breaking progress in numerical simulations over the 
last decade or so [TT1 - H4"] . More or less cuspy profiles towards the galactic centre result 
in significantly different J— factors, so it has become customary to use NFW, Einasto and 
isothermal profiles in order to bracket dark matter constraints. This sort of uncertainty, 
together with the local dark matter density that is uncertain to tens of percent [3"oT44"2"] . 
translates directly into the normalisation of the fluxes in figure [TJ Likewise, the particular 
choice of the observed target region affects sensitively the derived gamma-ray constraints 
(see for instance [251 H31 Si])- All this problematic has been exhaustively addressed in the 
literature; what we wish to motivate and pursue here is the search for a different kind of 
gamma-ray feature that is at least as powerful as the most well-studied spectral signatures. 
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III. RESULTS 



It is already evident from figure [TJthat thermal cross-sections for dark matter annihilations 
into intermediate scalars are in tension with the Fermi-LAT observations. We now turn to 
the derivation of precise constraints on the dark matter parameter space. A given dark 
matter model defined by (tudm, (o~v) , Am/moM) ~~ or ( m DM, I\ Am/mflM) - is considered 
excluded if the photon signal 7 plus background 7j h exceed the measurement more than 
2a in at least one Fermi-LAT energy bin over the range 1 — 280 GeV. In order to bracket 
the uncertainty regarding background modelling, we adopt three distinct approaches: (i) 
conservative, in which we boldly assume no background 7i b = 0, naturally leading to the 
weakest limits; (ii) intermediate, where we fit Fermi-LAT datapoints below 20 GeV to a 
power-law (scaled down to ~ 80% to give room to a putative dark matter signal) and take 
that to describe the gamma-ray background 7i & oc E~ u [y ~ 2.6); and (Hi) aggressive, 
where we take the background to meet the central points of the measurements which is a 
slightly overoptimistic method of setting constraints on dark matter models. Notice that the 
intermediate prescription represents an ad hoc approach of illustrative character and should 
be refined for more extensive analyses. This procedure, although not as sophisticated as 
others [23H25] , is appropriate for our purposes and enables an overall assessment of gamma- 
ray constraints on boxy spectra. 

Figure [2] (left) displays our constraints in the (av) vs m£>M plane for different values of 
Am/mpM and using our intermediate approach. The first overall impression is that the 
upper limits on the annihilation cross-section into a pair of scalars are rather strong, sitting 
up to two orders of magnitude below the usual thermal cross-section (av) = 3 x 10 -26 cm 3 /s. 
Considering our aggressive approach the constraints exclude even lower cross-sections of 
about 10~ 29 cm 3 /s at the lowest dark matter masses, as it is shown on the right panel 
of figure [2] This appears to be a general feature of cascade models and, as we shall see 
in Section |IV[ sets stringent limits on sensible particle physics models. The dips seen at 
itidm — 10 — 20 GeV in the left panel of figure [2] are artifacts of our intermediate approach 
and result from the energy bin at E n = 10 GeV where the measured flux is closest to our 
assumed background power-law. This makes the intermediate constraints for the almost 
degenerate case particularly strong at tudm — 20 GeV for annihilating dark matter (since 
E y = m DM /2) and m DM = 40 GeV for decaying dark matter (since E^ = m DM /^). The grey 
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FIG. 2: Gamma-ray constraints on cascade dark matter annihilations for the galactic centre 
region. In the left panel the lines show our constraints in the (av) vs moM plane for several values 
of the degeneracy parameter Am/rriDM and using the intermediate approach. In the right panel 
we compare the constraints obtained with the different aproaches for Am/moM = 0.001. The 
constraints obtained in Ref. [23] for 77 lines are represented by the filled circles. 



band in figure [2] (right) brackets the constraints derived with our conservative, intermediate 
and aggressive approaches. The band encompasses almost two orders of magnitude at low 
masses and less than one at high masses. This is because in the aggressive case the constraints 
are set by the Fermi-LAT error bars, which are much smaller at lower energies. Notice as 
well that, in the almost degenerate case at tudm ^ 600 GeV, the low-energy end of the box 
surpasses the last bin of Fermi-LAT data, leading to extremely weak constraints and to the 
sudden rise in the high-mass region of both panels in figure [2j 

Secondly, as expected, the strongest constraints are obtained for the almost degenerate 
case Am/mpM — 0.001, in which the gamma-ray spectrum resembles a line. In this case, 
a comparison with the results of Ref. |23j (shown by the filled circles in the right panel of 
figure [2]) is allowed, but the translation of the results is not immediate. In fact, a direct 
annihilation into 77 has a monochromatic spectrum at moM normalised to two photons per 
annihilation, while in the degenerate case of our model the line occurs at and is 

normalised to four photons per annihilation. In any case, as we can see on the right panel 
of figure [2j our aggressive approach gives similar results to the detailed likelihood method 
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presented in [23] . 

Finally, let us notice that, as Am/moM ~> 1> the constraints in figure [2] (left) weaken, 
but the weakening seems to "saturate" at the highest values of Amjm^M leading to rather 
robust constraints even in those cases. The reason behind such behaviour has already been 
addressed in Section [TTJ the decrease in the box amplitude as less and less degenerate cases are 
considered is compensated by the extension of the box edge E + to higher energies, where 
better sensitivities are achieved. This has the important consequence that, for virtually 
all relative degeneracies Am/mpM > 0.05, one can exclude roughly the same range of 
dark matter masses for fixed (cv), as more clearly seen in figure [3] (left). There we show 
the excluded regions of the dark matter parameter space for different values of (crv). For 
example, at (crv) = 3 x 10~ 26 (3 x 10~ 27 ) cm 3 /s and Am/mohi > 0.05 the dark matter 
mass range tudm — 5 — 600 (5 — 50) GeV is excluded using our intermediate approach. For 
completeness we show in figure [3] (right) the excluded parameter space in the case of decays 
cf) jZ, where the upper left white region is forbidden for kinematical reasons. In practice, 
even for particle physics realisations featuring highly non-degenerate configurations, it is 
possible to derive meaningful, stringent constraints. 

Figure [4] displays our intermediate constraints on the decay timescale for decaying 
dark matter taking Fermi-LAT observations of the halo region. As in the case of annihi- 
lations, it is clear that cascade decays into photons place relatively strong limits on dark 
matter models. In particular, the lower limits on T _1 fall in the range 10 28 — 10 29 s , i.e. two 
to three orders of magnitude above the typical 10 26 s used in the literature [45J to meet the 
PAMELA positron excess. The bumps at ttidm = 20 — 40 GeV present in figure [4] have the 
same origin as the dips in figure [2] (left) which have been explained above. 

Up to this point we have been focussing on the constraints provided by the search for 
box-shaped gamma-ray features in Fermi-LAT data. Let us now briefly comment on the 
detectability of these type of signatures. In the case of the widely-studied gamma-ray lines, 
a positive, statistically significant detection would give the peak energy and amplitude of 
the line which would translate directly into ttldm and (erf) (modulo background modelling 
and uncertainties regarding the J— factor). In the case of a box spectrum, instead, one 
would probably just spot the high-energy shoulder and accordingly estimate the box edge 
E + and amplitude. These two observables are obviously degenerate in the parameter space 
( m DM, ((^v), Am/m DM ). Thus, the discovery of a hard shoulder would be somewhat prob- 
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FIG. 3: Annihilating dark matter constraints in the Am/mpM vs tudm plane for different values 
of (av) using the intermediate approach. The shaded areas are excluded for the respective cross- 
sections. In the left (right) panel we show the case of intermediary particles decaying into 77 

m. 

lematic in reconstructing the corresponding dark matter phenomenological parameters and 
in connecting to concrete particle physics models. Nevertheless, it rests as a possibility that 
such spectral feature is detected at very high energies. 



IV. CONCRETE MODELS 



We consider a model where the dark matter is constituted by complex scalar particles, 
X, which are singlets under the Standard Model gauge group. To ensure the stability of the 
dark matter particle we introduce a discrete Z 2 symmetry, under which x is °dd while all 
the Standard Model particles are even. The Lagrangian then reads C = Csm + £>x + ^intj 
where C x contains the kinetic term and the potential term of the extra scalar x an d Ant is 
the interaction term of x an d the Standard Model particles, which simply reads: 

-£ int = \\X?\H\ 2 , (5) 

where H is the Standard Model Higgs doublet, H = ~^(0,v + h), with v = 246 GeV. In 
this simple model, the dark matter particle can annihilate XX ~^ hh with a sizable cross- 
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FIG. 4: Gamma-ray constraints on cascade dark matter decays for the halo region. The lines 
show our constraints in the T _1 vs tubm plane for several values of the degeneracy parameter 
Am/niDM and using the intermediate approach. 

section [46J. Subsequently, the Higgs particle h can decay h — > 77, with a branching ratio 
which is ~ 2x for the mass window of the light Higgs ~ 110 — 130 GeV [17], presently 
allowed by searches at the Large Hadron Collider [HI H9j. In this case, the branching ratio 
h — > 77 is too small to allow the observation of a gamma-ray signal, if the annihilation 
cross-section coincides with the thermal cross-section (av) ~ 3 x 10 _26 cm 3 /s. Nevertheless, 
it follows from figure [2] that a signal could be observed if the boost factor for XX ~~ >" hh 
is (9(20). Boost factors of this order are not unfeasible if a combination of substructure 
[50] . Sommerfeld enhancement [51] and/or non-standard cosmological scenarios [52] is at 
play. Furthermore, this limit on the boost factor is similar to the limit obtained in [53] 
for annihilations into Higgses from the first year of Fermi-LAT gamma-ray observations in 
the full sky, although it is likely weaker when compared to dwarf [51] and halo [55] recent 
observations. This simple example shows the power of the search for box-shaped gamma-ray 
features in constraining dark matter models. 

The intensity of the gamma-ray box signal can be enhanced by considering the following 
extension of the previous model. In addition to the dark matter particle \i we introduce a 
second, real, scalar field 0, singlet under the Standard Model gauge group and two fermions 
£1 and £ 2 , singlets under SU(3) C x SU{2) L and with hypercharge Y — 1 and —1 respectively 
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(the model is then free of gauge anomalies). Furthermore, we assume that <fi is even under 
the discrete Z 2 symmetry, while £1 and £2 are odd. 2 The Lagrangian then reads C = 
£sm + £ x + £</, + £g + Ant- Here, C x and are the Lagrangians for a complex and real 
scalar fields, respectively, contains the kinetic terms for the fermion fields £1 and £2, and 
dot contains the interaction terms among fields: 

-C mt = y£i£ 2 + + e£i(A M - tanOwZ^^ - e£ 2 (^ M - tan^Z^^ + h.c. 

+ ^A 0// 2 |#f + ^\ x h\x\ 2 \H\ 2 + ^A 0x 2 |x| 2 + ^\H\ 2 + «0|x| 2 • 

With this interaction Lagrangian, the dark matter particle \ can annihilate into two scalar 
particles, XX ~~ >* 00> as wen as into hh, ZZ, WW or // [46J. Assuming that the coupling 
\$ x is much larger than X x h, then the dark matter particle predominantly annihilates into 
the intermediate 0-scalars. 

The annihilation is followed by the decay of (p. At tree level <fi can only decay into 
two Higgs particles <fi — > hh, provided the decay is kinematically allowed, namely when 
777,0 > 2m/j. In this case the decay rate is 



ro -> hh) 



128vr 



1 - 



Ami 



1/2 



(6) 



This is usually the dominant decay mode, unless the coupling /1 is very small. 

At the one loop level, on the other hand, new decay channels are possible and can have a 
sizable branching ratio when the decay — > hh is kinematically forbidden. Since £12 carry 
hypercharge, <fi can decay via an one loop diagram into 77 and, if kinematically allowed, 
also into 7Z and ZZ. The rates for the decay processes into gauge bosons read, in the limit 
m € > M z , rrifa 



r(0 -»■ 77) 
r(0 1 z) 

r(0 -> zz) 



1447T 3 m? 



y 2 a: 2 tan 2 W m 



72tt 3 



y 2 a 2 tan 4 6» w m 

1447T 3 




1 - 



Ml 



mz 



4M| 

2 

mz 



1/2 



2M| 

2 

mi 



+ 



2M| 

2 

mz 



(7) 



2 Assigning an odd charge to under Z 2 is not determined by the dark matter physics, but to prevent the 
mixing of the right-handed electron with £ 2 - 
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Furthermore, the scalar can decay — > bb via a loop diagram involving the coupling 
^ficf)hh to the Higgs particle. The decay rate for this process reads 



r(0 66) 
where, taking m& <C m^, 

/M = 



128vr 3 V m 2 



4m2 



1 /-l-a; 

cfe / oh/ 
o Jo 



ft / 

1 — x — y 



(8) 



(9) 



x + y — Awxy 

We show in figure [5] the branching ratios for decays in this model, assuming for con- 
creteness y — 1, = 500 GeV, /x = 100 GeV and = 125 GeV. As apparent from the 
figure, when < 2m/,, the scalar decays dominantly into two photons. Therefore, in this 
situation the stringent constraints in the parameter space presented in figure [3] (left) apply. 
Let us note however that in this particular model the dark matter particle is not Majorana, 
and accordingly the flux in equation ^ is to be multiplied by 1/2. 
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FIG. 5: Branching ratios for the decay of the intermediate scalar particle for the concrete model 
described in the text. 

The study of more complex particle physics models in which box constraints are relevant 
is deferred to future work. 



V. CONCLUSION 



In this work we have paid close attention to a class of gamma-ray spectral signatures that, 
to the best of our knowledge, has not been explored in the framework of dark matter searches. 
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Cascade dark matter annihilations or decays into photons give rise to a box-shaped gamma- 
ray spectrum that resembles a line in the limit of degenerate dark matter and intermediate 
scalar masses, and a flat plateau in the non-degenerate case. This is the feature of interest 
for the present paper. The key point is that, as the mass splitting increases, the photon 
spectrum gets dimmer but reaches higher energies, yielding sizable fluxes even for highly 
non-degenerate configurations. We have illustrated how this unique phenomenology applies 
to dark matter searches by using the 2-yr Fermi-LAT data. Overall, our constraints for 
both annihilating and decaying dark matter are rather strong, reaching (av) ~ 1CT 29 cm 3 /s 
(i.e. ~ 3 x 10" 4 in branching ratio assuming a thermal cross-section) and T" 1 ~ 3 x 10 29 s 
in the most aggressive cases. Also, we have found a rapid "saturation" of the constraints 
as the mass splitting between dark matter and intermediate scalar is increased. Concretely, 
any mass splitting Am/moM above ~5% leads essentially to the same annihilation cross- 
section upper limits, which are only a few times weaker than the degenerate configuration. 
A similar situation holds for the decaying dark matter constraints. This means that no 
fine-tuning between the masses at play is required to place stringent limits on the properties 
of dark matter. Consequently, the non-observation of box-shaped gamma-ray features is 
extremely effective in constraining concrete particle physics frameworks. We point out that 
more detailed research regarding this alternative spectral signature is needed in view of the 
exquisite data provided by the Fermi-LAT instrument as well as the higher energy data from 
IACTs and, eventually, CTA. 
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